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The electronic structures and optical properties of cubic HfOg are calculated by gen- 
eralized gradient approximation (GGA)+[/ approach. Without on-site Coulomb in- 
teractions, the band gap of cubic HfOj is 2.92 eV, much lower than the experimental 
value (5.7 eV). Introducing the Coulomb interactions of 5d orbitals on Hf atom ([/'') 
and of 2p orbitals on O atom {U^), we can reproduce the experimental value of the 
band gap. The calculated dielectric function of cubic Hf02 by the GGK+V^+U^ 
approach predicts the presence of a shoulder structure below the main peak of the 
absorption spectrum. These indicate that the QQK+U'^+U^ approach is a conve- 
nient and powerful method to calculate and predict the electronic structures and 
optical properties of wide-gap optical materials. 
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I. INTRODUCTION 



Hf02 has been widely studied both experimentally and theoretically because of its excel- 
lent dielectric properties, wide band gap, high melting point, etc.-i^. Moreover, HfOg has 
been proved to be one of the most promising high-dielectric-constant materials, due to its 
high formation heat and nice thermodynamic stability when it contacts indirectly with Si^. 
It has also been used in optical^ and protective coatings^. At atmospheric pressure, Hf02 
exists in three polymorphs when temperature changes^: monoclinic structure (P21/c, C2h 
5) at low temperatures, tetragonal (P42/nmc, D4h 15) at temperature higher than 2000 K, 
and finally, with increasing temperature over 2870 K, the cubic fluorite (Fm3m, Oh 5) be- 
comes stable. However, by rapid quenching or stabilizer addition, the cubic structure can 
be stabilized at low temperatures, e.g., see Refs. Q and 8. 

The electronic and structural properties of HfOg have been theoretically studied by the 
first-principles band structure calculation based on the density-functional theory (DFT). 
For example, the structure, vibration and lattice dielectric properties of Hf02 have been 
investigated in the frame of the local density approximation (LDA) and generalized gradient 
approximation (GGA)^. Based on the ABINIT package, the Born effective charge tensors, 
phonon frequencies and dielectric permittivity tensors in cubic and tetragonal phases are 
calculated by LDA—. The full potential linearized augmented plane-wave method based on 
the WIEN2k code is also employed to study the structural and electronic properties in the 
cubic phased. Further, the spin-orbit effect is estimated within the framework of LDA and 
GGA— . The resulting band gap is about 3.65 eV, still smaller than the experimental value 
(5.7 eV)^, indicating that correlation effects in HfOg cannot be ignored. 

Jiang et alr^ have studied the band structures of Hf02 by the GW approximation, where 
the self-energy of a many-body system is approximated in terms of the screened Coulomb 
interaction (W) up to first order, and then correlation effects can be addressed to some 
extent. The resulting band gap value is 5.2 eV. Although the value is much closer to the 
experimental one, the GW method is highly cost in terms of numerical demands. Another 
technique to include correlation effects with relatively less computational effort is the so- 
called LDA+f/ or GGA+f/ approach, in which the correlation effect is incorporated through 
the on-site Coulomb interaction U—. As an approximation to GW, GGA(LDA)+t/ approach 
can reproduce experimental data as accurately as the hybrid functional approach but with 



much lower computational efforts^. It can give a qualitative improvement compared with 
the LDA, not only for excited-state properties such as energy gaps, but also for ground-state 
properties such as magnetic moments and interatomic exchange parameters'^. Recently, a 
GGA+f/ approach for CeOg has been performed, where not only the on-site U for Ce 4/ 
electrons, but also for O 2p electrons, are included^. The choice of U is not unambiguous. 
Though there are attempts to extract it from the first-principles calculations, it is nontrivial 
to determine its value a priori. Hence, in practice, U is often fitted to reproduce a certain 
set of experimental data, such as band gaps and structural properties. 

In this paper, we use GGA and GGA+U schemes formulated by Loschen et al.—, to 
calculate the lattice parameters, band structures, and optical properties of cubic HfOg (c- 
Hf02). We find that the value of the band gap can be reproduced by introducing the on-site 
Coulomb interactions of 5d orbitals on Hf atom {W^) and of 2p orbitals on O atom (U^). 
We also notice that a shoulder structure below the main peak of the imaginary part of 
the dielectric constant remains even when and are introduced. This implies that 
the shoulder structure is robust against the change of band structure due to the Coulomb 
interactions. We thus expect that the shoulder will appear in absorption measurements for 
pure c-IIf02. These findings imply that the GGA+W^+U^ approach is a convenient and 
powerful method to calculate and predict the electronic structures and optical properties of 
wide-gap optical materials. 

II. COMPUTATIONAL METHODOLOGY 

Density functional theory calculations are performed with plane-wave ultrasoft pseu- 
dopotential, by using GGA with Perdew-Burke-Ernzerhof (PBE) functional and GGA-|-?7 
approach as implemented in the CASTEP code (Cambridge Sequential Total Energy Pack- 
age)^. The ionic cores are represented by ultrasoft pseudopotentials for Hf and O atoms. For 
Hf atom, the configuration is [Xe]4/^^5(i^6s^, where the 5(P and 6s^ electrons are explicitly 
treated as valence electrons. For O atom, the configuration is [He]2s^2p^, and valence elec- 
trons include 2s^ and 2p^. The plane-wave cut off energy is 380 eV. And the Brillouin-zone 
integration is performed over the 24 x 24 x 24 grid sizes using the Monkorst-Pack method 
for cubic structure optimization. This set of parameters assure the total energy convergence 
of 5.0 X lO^^eV/atom, the maximum force of 0.01 eV/A, the maximum stress of 0.02 GPa 
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and the maximum displacement of 5.0 x 10 ^A. 

In the following sections, we firstly optimize the geometry structure of c-HfOg by the 
GGA method. Using the optimized structure, we next introduce If^ for Hf 5d orbitals and 

for O 2p orbitals. Comparing the numerical values of the band gap with the experimental 
one, we obtain the best values of and U^. The electronic structures and optical proper- 
ties of c-HfOg are calculated by means of GGA, without U and with W^+U^, respectively. 
Comparison with available experimental data is presented. 



III. RESULTS AND DISCUSSION 
A. GGA Calculation 

The space group of c-HfOg is Fm3m and the local symmetry is 05h. Moreover, c-HfOg is 
fully characterized by a single lattice constant a. The GGA calculation of the perfect bulk 
c-HfOg is performed to determine an optimized a in order to check the applicability and 
accuracy of the ultrasoft pseudopotential. The optimized a is 0.526 nm, in good agreement 
with other theoretical values^"— The deviation from experimental values (0.508 nnt^S and 
0.516 nn>i^) is 3.54% at most. Hereafter, we use this optimized value of a. 

The band structure along high-symmetry directions of the Brillouin zone and total density 



of states (DOS) of c-Hf02 are shown in Fig. [TJ The band structure in Fig. 1(a) shows a 
direct band gap because the top of the valence bands and the bottom of the conduction 
bands are found at the same X point. The value of the band gap Eg is around 2.92 eV, much 
smaller than the experimental value 5.7 eV. This is due to the well-known underestimate 
of conduction-band energies in ab initio calculations: the DFT results often undervalue the 
energy of 5d orbitals of Hf atom, lowering the bottom level of conduction bands. As a result, 
the band gap obtained by GGA is lower than the experimental one. 



The total DOS is presented in Fig. 1(b) There are two parts in the valence band, namely. 



the lower region from — 18.9eV to — 15.2eV and the upper region from — 6.6eV to 0.1 eV. 



Figure 1(d) shows that the lower valence bands are predominantly composed of O 2s, while 



the upper valence bands consist of O 2p accompanied with hybridization with Hf 5d as 



shown in Fig. 1(c) The conduction bands below 10 eV are mostly composed of Hf 5d with 



some amount of O 2p. The 6s and 5p orbitals of Hf atom also contribute to the conduction 



-20 J 1 1 

W L G X W K 



-20 C— 1 ■ 1 ■ 1 ■ 

1 2 3 4 5 

Total DOS 



(a) 



(b) 



CO 

O 
a 



Hf 




-10 10 20 

Energy (eV) 



m 2 - 
O 

Q 



-10 10 20 

Energy (eV) 



(c) (d) 

FIG. 1. The band structure and density of states (DOS) of c-HfOg obtained by GGA. (a) Band 
structure; (b) total DOS; (c) partial DOS of Hf atom, and (d) partial DOS of O atom. 



bands, though their values are small compared with 5d states. 

Dielectric function is obtained by taking into account of inter-band transitions. Figure [2] 
shows the complex dielectric function as a function of photon energy. The real part of the 
dielectric function ei in the low energy increases with photon energy and gets maximum at 
3.5 eV. Then it drops sharply with photon energy to negative values and return back to 
positive later. The imaginary part £2 shows a maximum at 7.5 eV followed by two shoulder 
structures below the energy. The two shoulders may be related to DOS in conduction bands 



at 4eV and 7eV as shown in Fig. 1(a) Our results of the dielectric function agree with 
other calculations for c-HfO^"^. However, the magnitude of maximum value of £2 (~ 10) 
is larger than an available experimental value (~ 8) of HfO^. This disagreement will be 
resolved if the gap magnitude is reproduced by introducing the on-site Coulomb interactions, 
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FIG. 3. Calculated band gap Eg as a function of (a) C/<^ and (b) UP. 
as discussed below. 



B. GGA+U'^+UP Calculation 

Using the optimized lattice parameter, a = 0.526 nm, we calculate the band structure 
and DOS of c-Hf02 by changing If^ for 5d orbitals of Hf atom. The band gap Eg obtained 



from the band structure is shown in Fig. 3( function of U . It can be seen that E, 
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firstly increases, and then drops with increasing U , showing a maximum value (4.43 eV) 



at = 8.25 eV. The maximum value is smaller than the experimental one (5.7 eV). The 
saturation of Eg with may be related to the approach of 5d states toward 6s and 5p 
states, though microscopic mechanism is not yet fully understood. Next, we introduce 



for 2p orbital of O atom, while keeping U = 8.25 eV. The result in Fig. 3(b) shows that 



Eg monotonically increases with UP. When U'^ = 8.25 eV and = 6.25 eV, the calculated 
band gap is 5.697 eV, coinciding with the experiment one. 

There could of course exist different combinations of and which can reproduce the 
band gap. For example, = 10 eV, = 6.35 eV, the resulting band gap is 5.702 eV, with 
similar optical properties in the low-energy regime. Based on physical consideration with 
reference to other numerical results, here, we choose If^ = 8.25 eV, = 6.25 eV as a typical 
representative and perform the GGA+[/ calculations. 



The resulting band structure is shown in Fig. 4(a) The bottom of the conduction bands 



is moved to the G point. The band dispersions near the bottom shift to higher energy with 



the increase of If^ and are reconstructed as compared with the dispersions in Fig. 1(a) As 



a result, the separated DOS at 4eV and 7eV in Fig. 1(b) merges to one sharp structure at 



6eV in Fig. 4(b) It is clear from Fig. 4(c) that the reconstruction is caused by the 5d state 



of Hf. The O 2p states are also affected by the reconstruction through hybridization effect 



as seen in Fig. 4(d) We note that the band gap increases with increasing through such 
strong hybridization effect. 

Figure [5] shows the dielectric function. The real part ei exhibits a maximum at 6.38 eV, 
corresponding to band edge reflection peak at 5.0 eV in the reflection spectrum. The calcu- 
lated static dielectric constant is 3.44, coinciding with the experimental value of monoclinic 
HfO^i^. The imaginary part €2 shows a maximum at 8 eV. The maximum value (~ 7.6) 
is very close to the experimental value (~ 8)^, in contrast with the GGA case as men- 
tioned before. This improvement comes from the agreement of Eg between the numerical 
and experimental values. Another interesting observation in €2 is the presence of a shoulder 
structure below the peak. It is remarkable that, even though the reconstruction of the band 
structure occurs near the bottom of the conduction band, the shoulder structure remains as 
is the case of the GGA calculation. This implies the robustness of the shoulder structure in 
c-Hf02. Therefore, we can expect that such a shoulder structure appears in the absorption 
spectrum in pure c-Hf02 samples. 

Other optical properties can be computed from the complex dielectric function^^. We have 
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FIG. 4. The band structure and density of states (DOS) of c-Hf02 obtained by GGA + If^ 
{V^ = 8.25 eV) + UP {U'p = 6.25 eV). (a) Band structure; (b) total DOS; (c) partial DOS of Hf 
atom, and (d) partial DOS of O atom. 



calculated the refractive index of c-HfOg by GGA+f/'^+f/^, including refractive coefficient n 
and extinction coefficient k. The results show that n is 1.85, consistent with the experimental 



value of monoclinic HfOg, 1.93^. 



IV. CONCLUSIONS 

The structure optimization of c-Hf02 is performed by using first-principles GGA and the 
resulting cell parameter of optimization structure is 0.526 nm, consistent with the experi- 
mental value and other theoretical results. However, the numerical value of band gap is 
only 2.92 eV, much lower than the experimental ones. The deviation can be resolved by 
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FIG. 5. Real and imaginary parts of dielectric function for c-Hf02 obtained by GGA+W^+U^. 

introducing on-site Coulomb interactions into 5d orbital of Hf atom (^7'^) and 2p orbital of 
O atom ([/^) at the same time. We obtain the best estimations of and for c-Hf02. 
The electron structure and optical properties are calculated both by GGA and GGA + If^ 
{W^ = 8.25 eV)+ JJP {U^ = 6.25 eV). We find that a shoulder structure below the main 
peak of the imaginary part of the dielectric constant remains even when and are 
introduced. This implies that the shoulder structure is robust, and we thus expect that 
the shoulder appears in the optical measurements for pure c-Hf02 samples. These findings 
imply that the GGA+W^+U^ approach is a powerful method to calculate and predict the 
electronic structures and optical properties of wide-gap optical materials. 
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